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ABSTRACT 

Context. VISTA variables in the Via Lactea is an ESO Public survey dedicated to scan the bulge and an adjacent portion of the 
Galactic disk in the fourth quadrant using the VISTA telescope and the near-infrared camera VIRCAM. One of the leading goals of 
the VVV survey is to contribute to the knowledge of the star cluster population of the Milky Way. 

Aims. To improve the census of the Galactic star clusters, we performed a systematic and careful scan of the JHK^ images of the 
Galactic plane section of the VVV survey. 

Methods. Our detection procedure is based on a combination of superficial density maps and visual inspection of promising features 
in the H-, and Kg-band images. The material examined are VVV JHKs color-composite images corresponding to the Data 
Release 1 of VVV. 

Results. We report the discovery of 493 new infrared star cluster candidates. The analysis of the spatial distribution show that 
the clusters are very concentrated in the Galactic plane, presenting some local maxima around the position of large star-forming 
complexes, such as G305, RCW 95, and RCW 106. The vast majority of the new star cluster candidates are quite compact and 
generally surrounded by bright and/or dark nebulosities. IRAS point sources are associated with 59% of the sample, while 88% are 
associated with MSX point sources. GLIMPSE 8 pm images of the cluster candidates show a variety of morphologies, with 292 
clusters dominated by knotty sources, while 361 clusters show some kind of nebulosity in this wavelength regime. Spatial cross¬ 
correlation with young stellar objects, masers, and extended green-object catalogs suggest that a large sample of the new cluster 
candidates are extremely young. In particular, 104 star clusters associated to methanol masers are excellent candidates for ongoing 
massive star formation. Also, there is a special set of sixteen cluster candidates that present clear signspot of star-forming activity 
having associated simultaneosly dark nebulae, young stellar objects, extended green objects, and masers. 

Key words, open clusters and associations - galaxy: stellar content - galaxy: structure - infrared: stars - surveys 


1. Introduction 

Stellar clusters play a fundamental role in the dynamics and 
chemical evolution of galaxies. They trace recent star formation 
history and display some of the most spectacular scenarios of 
stellar evolution. Indeed, knowing how the clustered stellar pop¬ 
ulations are distributed within our galaxy is essential in the study 
of its large scale structure, hence a key objective of large sur¬ 
veys like VISTA variables in the Via Lactea (VVV, Minniti et 
al. 2010). However, newborn stars formed in clustered environ¬ 
ments are often embedded in large or giant molecular clouds that 


* Based on observations taken within the ESO VISTA Public Sur¬ 
vey VVV, Program ID 179.B-2002. The full catalog is only available at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsarc.u-strasbg.fr/viz-bin/qcat7J/A-hA/. 
rbarba@userena.cl 


generate heavy obscuration, which hinder optical observations of 
their early development stages, so observations at infrared wave¬ 
lengths are preferred (Lada & Lada 2003). 

In the particular case of the Milky Way, any study of its 
global structure and morphology is difficult because we are look¬ 
ing at it from a location placed just in the Galactic plane, where 
the bulk of the gas, dust, and stars are known to be confined, 
making the extinction and crowding extremely complex and 
high. Because of that, and in spite of developments in instrumen¬ 
tation and observational techniques, a significant fraction of the 
Milk Way cluster population remains unknown. In fact, it has 
been estimated that our galaxy should have at least ~ 23,000^ 


^ This number is a lower limit for the estimated population of star clus¬ 
ters in the Galaxy. Bonatto et al. (2006) indicated a total number of open 
clusters in the range 1.8-3.7x10^. 
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open clusters (Portegies Zwart et al. 2010), however only a few 
thousand of them have been identified to date. 

Customarily, the detections of new Galactic star clusters are 
mostly based on a visual inspection of recorded photographic 
and electronic images (CCDs, IR detectors, etc.) taken both at 
optical and near-infrared (NIR) wavelengths, and the develop¬ 
ment of new homogeneous surveys in the past decade has al¬ 
lowed the discovery and analysis of many unknown cluster can¬ 
didates: e.g., Bica et al. (2003), Dutra et al. (2003), Mercer 
et al. (2005), Froebrich, Scholz, & Raftery (2007), Borissova 
et al. (2011), Solin et al. (2014), etc. Particularly relevant to 
our research is the Two-Micron All-Sky Survey (2MASS, Skrut- 
skie et al. 2006), the prior near-infrared all sky survey, covering 
the same wavelength range as VVV, but different spatial resolu¬ 
tion (4'.'0 vs.O'.'?) and photometric depth (^s-mag 15.5 vs. 18 in 
VVV). Also, if we consider that open clusters are gravitationally 
bound structures that consist of more or less coeval stellar popu¬ 
lations, the detection of density enhancements, either visually or 
by means of computational algorithms, should only be the first 
step to properly identifying potential new clusters. Indeed, stel¬ 
lar density enhancements may also be produced by chance align¬ 
ments or localized low foreground extinction (e.g., Odenkirchen 
& Soubiran 2002; Maciejewski & Niedzielski 2008; Moni Bidin 
et al. 2010; Froebrich, Scholz, & Raftery 2007). In this sense, to 
confirm that a group of stars are really physically related and 
forming a true cluster, some additional criteria such as radial 
density profiles (e.g., Gaussian or King), colors, and/or kinemat¬ 
ics, are probably necessary. Besides that, the Hipparcos catalog 
(ESA 1997) has provided accurate positions (as probably GAIA 
will do in the future), proper motions, and trigonometric paral¬ 
laxes for a significant stellar population in the solar neighbor¬ 
hood. For example, many OB associations show small velocity 
dispersions (e.g., Tian et al. 1996), which makes them appear 
as coherent structures in the velocity space, and de Zeeuw et al. 
(1999) used this characteristic to distinguish the membership in 
the Hipparcos catalog of nearby OB associations inside a ~ 1 kpc 
radius of the solar neighborhood. 

Here we present the results of our search for new cluster can¬ 
didates in the Galactic disk portion covered by the VVV Data 
Release 1 (DRl) images (Saito et al. 2012). This work is based 
on superficial density maps and visual inspection of the associ¬ 
ated H-, and Ks-band VVV DRl images. As a result almost 
500 new Galactic cluster candidates were identified in the disk 
area of the VVV survey. To make this catalog, we performed a 
more systematic search in each VVV tile than Borissova et al. 
(2011) using the same set of VVV DRl images. This paper is 
organized as follow: In Section 2, we describe briefly the main 
characteristics of the VVV. In Section 3, we describe the find¬ 
ing methods for the new star cluster candidates. In Section 4, 
we describe the characteristics of the new star clusters and the 
spatial cross-correlation with different infrared and radio source 
catalogs, and in the Section 5, we summarize our conclusions. 

2. Data: VVV tile images 

The VVV survey observations are being carried out in the four- 
meter Visible and Infrared Survey Telescope for Astronomy 
(VISTA), using the VIRCAM (VISTA InfraRed CAMera; Dal¬ 
ton et al. 2006; Emerson & Sutherland 2010) located at ESO’s 
Cerro Paranal Observatory in Chile. VIRCAM has a field of view 
of 1.65° by side, but the mosaic of 4 x 4 Raytheon 2048 x 2048 
pixels^ detectors only covers 0.6 deg^ in each single exposure (a 
pawprint), with spacings of 40% and 90% in the X and Y axes, 
respectively, and a pixel scale of O'.'34. Typically, six VISTA 


pointings (pawprints) are combined to form a complete VVV 
field, a tile. The VVV survey was planned to observe 348 of 
these fields in the Galactic bulge and a portion of the Southern 
Galactic disk (196 and 152 tiles, respectively), in five broadband 
filters Z, F, /, //, Ks. The area covered by the VVV survey over¬ 
laps completely with the 2MASS survey, with part of the Spitzer 
Space Telescope GLIMPSE survey (Benjamin et al. 2003), and 
also complements the area covered by the UKIDSS survey (Lu¬ 
cas et al. 2008). 

This work is based on the analysis of the 152 /, //, and Ks 
Galactic plane tile images, which were produced by the Cam¬ 
bridge Astronomical Survey Unit (CASU) from the scheduled 
observations of the first and second years of the survey. These 
images are available as DRl (Saito et al. 2012). Further details 
about the VVV survey, technical description, observing strategy, 
and scientific goals are in Minniti et al. (2010). 

3. Star cluster candidate finding methods 

Our approach for detecting and selecting each new star cluster 
candidate was to look for variations in the stellar density by 
simple star counts, helped by the systematic examination of the 
tile image of each band and the combined color mosaics dis¬ 
played in the Aladin virtual observatory tool (Bonnarel et al. 
2000). The selected color combination was the usual: /, H, and 
Ks images combined as blue, green, and red channels, respec¬ 
tively. The star-count methodology is well known and has been 
used by several researchers (e.g., Lada & Lada 1995; Carpenter, 
Snell, & Schloerb 1995; Carpenter, Heyer, & Snell 2000; Ku¬ 
mar, Kamath, & Davies 2004; Kumar, Keto, & Clerkin 2006). 
The method can be refined by using different bin sizes in order 
to investigate large-scale structures, as well as smaller scale sub¬ 
clustering (e.g., Kumar, Kamath, & Davies 2004; Kirsanova et 
al. 2008), and/or also by smoothing the binned data over adja¬ 
cent bins (e.g., Karampelas et al. 2009). 

The star-count technique was applied to the Galactic disk 
portion of the VVV DRl Ks-band images. The stars were de¬ 
tected after using of the task DAOFIND (Package Daophot) in 
IRAF^. The limiting magnitudes are those for point sources pre¬ 
senting minimum peak count values, of 3cr above the local sky, 
which in general correspond to /, //, and Ks magnitude values 
of about 19.5, 18.8, and 18, respectively. Of course, such values 
must be considered as mean values because the numbers depend 
critically on the position of the sources relative to zones of ex¬ 
tended emission, as well as to the heavily crowded zones like 
those in the Galactic plane. It consists in dividing each tile into 
small bins of equal size and determining the number of stars in 
each bin. As result, those bins presenting values greater than 3cr 
above what was measured in a given control region were con¬ 
sidered as locations of potential candidate clusters, deserving 
further “by eye” examination. The binning sizes were chosen 
carefully such that the number of objects per bin was neither too 
small (prohibiting a meaningful analysis) nor too large (hiding 
existing features). 

We used tipical bin sizes varying from 30 pixels (about 10") 
to about 90 pixels (about 30"), which were determined empiri¬ 
cally on known cluster candidates in the Hit regions RCW 106 
and RCW 95, which were previously studied in the past by one 
of us (ARL). 

^ IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 
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Fig. 1: RGB color images of a representative sample of the newly discovered cluster candidates. The field-of-view for the chart of 
each cluster is set accordingly between 1x1 aremin, 3x3 aremin, and 5 x 5 in order to visualize the cluster comfortably. This set 
of clusters are the same showed in Figure 9. The complete set of charts for the whole catalog is offered in electronic form. 


As a quantitative example of the procedure applied in the de¬ 
tection of a local maximum in the superficial density maps on 
the ^s-band VVV images, for a “small” cluster candidate like 
La Serena 004 (with an estimated radius smaller than 30"), the 
binning size that enabled us to detect it with a appropriate con¬ 
trast against the surrounding field stars was 5". On the other 
hand, for a “large” cluster candidate like La Serena 012, the 
best results were achieved when we applied the biggest binning 


size of 30". Of course, we must point out that using this binning 
range possibly led to a cluster candidate detection rate biased 
toward small-to-medium sized cluster candidates, with most of 
the larger and sparse clusters being hard to detect. The control 
regions used in determining of the mean background to be sub¬ 
tracted from the statistics for each cluster candidate were com¬ 
puted from concentric rings with a mean radius equivalent to two 
to three times the observed estimated size of the respective clus- 


Article number, page 3 of 13 




























A&A proofs: manuscript no. Almost-500-cluster-candidates-VVV 


ter candidate. We also tested the use of adjacent control fields 
placed at several arc minutes (typically 2-5) from the center of 
the cluster candidate, finding no significative differences from 
each method. 

Finally, surface density maps were then compared with the 
color-composite JHKs images, which in turn were used to prop¬ 
erly evaluate and identify (by visual inspection of the local con¬ 
centrations of stars) those features related to candidate clusters. 
Many cluster candidates have affected the star counts owing to 
the presence of saturated stars. We must take into account that in 
the VVV images, the saturation limit in the ^s-band images is 
about magnitude 12.5. Additionally, we visually identified some 
compact stellar groups mostly associated with bright and/or dark 
nebulosities. If such a close association was confirmed, the clus¬ 
ter candidate was evaluated for overdensities and included in 
the catalog. Also, during the search, we discovered some con¬ 
centrations of saturated bright sources in the VVV images, but 
their appearance resembles a cluster of red giant or supergiant 
stars, similar in morphology to some clusters discovered in the 
2MASS survey (cf. Figer et al. 2006; Clark, J. S. et al. 2009). 
We included some of these groups after considering the 2MASS 
data for the saturated stars, for example, the star cluster candi¬ 
date labeled as La Serena 109. We checked carefully that each 
stellar group candidate had been not discovered previously, and 
any previously known neighbor cluster should be located as at 
the minimum distance of one arcminute. 

Catalogs of star clusters and star cluster candidates used for 
the check are Bica et al. (2003), Dutra et al. (2003), Mercer et 
al. (2005), Froebrich, Scholz, & Raftery (2007), and Borissova 
et al. (2011). Recently, Solin et al. (2014) have published a new 
catalog with star cluster candidates, using the same set of the 
Galactic disk VVV DRl images and also including the bulge 
area. These authors report 88 new star cluster candidates and 
39 star formation location candidates detected by a fit of a mix¬ 
ture model of Gaussian densities and background noise, and the 
expectation maximization algorithm to prefiltered near-infrared 
survey stellar catalog data. 

Interestingly, this independent search of cluster candidates is 
a very good alternative to checking the reliability of our meth¬ 
ods. In the VVV tiles of the Galactic disk, Solin et al. (2014) de¬ 
tect 82 star cluster candidates, 34 star formation location candi¬ 
dates, and 24 infrared nebulae of unknown nature. We performed 
a cross-match with a radius of 60" between our catalog and Solin 
et al. (2014) sources, obtaining a very good agreement. We “re¬ 
cover” 111 sources from the complete list of 140 objects in the 
Solin et al. (2014), which are 105 sources inside the 15" cross¬ 
matching radius. The “recovered” objects from the Solin et al. 
(2014) correspond to 68 cluster candidates (83% of Solin et al. 
list), 25 star forming region candidates (73%), and 18 infrared 
nebulae (75%). 

In Table 1 we present a sample of the catalog, as an example 
of its form and content. The full catalog includes 493 star cluster 
candidates, and it will be hosted by CDS^ database and also the 
VVV-ULS project page. Characteristics listed in columns 7 to 
11 will be described in the following sections. The catalog is 
organized as follows: 

Column 1 lists the star cluster candidate designation La Serena 

and running numbers from 1 to 493. 

Columns 2 to 5 gives the equatorial (J2000.0 epoch) and galac¬ 
tic coordinates of the cluster center. These coordinates are 

estimated from the bin center, and in many cases they are 

corrected with inspection of the images by eye. 

^ http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 


Column 6 indicates the VVV tile where the cluster is located. 
In a few cases, when the cluster is close to the edge of the 
tile, it is possible to find it in two different tiles. 

Column 7 lists the IRAS source associated with the cluster. 
Column 8 lists the MSX source associated with the cluster. 
Column 9 indicates the morphologic appearance of the cluster: 
O (open cluster) or C (compact cluster). Compact clusters are 
small clumps of stars, with sizes less than 30", mostly associ¬ 
ated with compact nebulosities in the VVV and/or GLIMPSE 
8 pm images, while on the other hand, larger extended clus¬ 
ters are classified as open clusters. 

Column 10 describes the size of the clusters based on an upper- 
limit estimation of the number of sources (Vs) associated 
with the cluster. Star clusters with < 20 are classified 
as small (S). Those with 20 < Vs < 50 are classified as 
medium (M). Finally, clusters with Vs > 50 are classified as 
large (L). We estimate these upper limits from the star count 
in the area of the cluster, doing an eyeball estimation of the 
size of the cluster. These numbers are for guidance because 
it must be noted that no background or crowding correction 
were applied. 

Column 11 describes the morphology of each cluster based on 
8 pm band images obtained by the Spitzer Space Telescope, 
under GLIMPSE survey. There are four main morphological 
types: point sources (P), knots (K), bubbles (B or sB), and 
nebula (N). Additional comments about the morphology of 
the clusters are annotated in Table 2. Point sources refer di¬ 
rectly to a group of point sources without noticeable nebulos¬ 
ity. Knot refers to a remarkable and compact nebulosity asso¬ 
ciated with point sources in the 8 pm images. Bubble refers 
to arc-shaped nebulosity associated with the point sources. 
In some cases the arcs or filaments form small circles delim¬ 
iting small cavities in the dusty clouds, while in other cases, 
the arcs define larger bubbles. For the case of small bubbles, 
they are marked as sB. Nebula refers to diffuse nebula with¬ 
out a specific shape associated with the point sources. 

Also, we present in Table 2 a sample of the association of 
different astrophysical sources with each cluster candidate ob¬ 
tained in the cross-matching procedures described in the Section 
4.2. Again the full table will be available in the CDS database. 
Table 2 lists the Solin et al. (2014) object name (Column 2) and 
the different types of objects associated with the new star cluster 
candidates: dark clouds (Column 3), young stellar object candi¬ 
dates (Column 4), extended green objects and outflows (Column 
5), and masers, with the corresponding bibliographic references 
(Columns 6 and 7). 

Additionally, in the Table 3 are included special notes and 
comments for each cluster candidate. 

We have prepared a website with the detailed astromet¬ 
ric and morphologic information for the 493 cluster candi¬ 
dates, including color-composite charts. This catalog is available 
at http://astro.userena.cl/science/lsclusters.php 
and it is open to the community. These color-composite images 
were created by combining /-band, i/-band, and Ks-band im¬ 
ages following the usual scheme of blue, green, and red chan¬ 
nels. Each color image was cropped and dynamically stretched 
to display each cluster comfortably. For example, the chart for 
the compact cluster La Serena 309 was cropped to a size of 
1' X 1', while the chart corresponding to the open cluster can¬ 
didate La Serena 109 was cropped to a 3' x 3' field of view, re¬ 
spectively. As example. Figure 1 presents JHKs color-composite 
images of some of the newly discovered star cluster candidates. 
The complete set of color-composite chart images will be avail¬ 
able in VVV-ULS web page. 
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Fig. 2: Angular distribution of the newly discovered La Serena star cluster candidates in the Galactic plane (upper left panel). Also, 
the histograms of the distributions of clusters in Galactic latitude (upper right panel) and Galactic longitude (lower panel). 


4. Characteristics of the new ciuster candidates 

4.1. Spatial distribution 

The angular distribution of the new stellar group candidates (Fig¬ 
ure 2, upper left panel) shows an obvious concentration toward 
the Galactic plane. The distribution in Galactic latitude is char¬ 
acterized by a single-peaked distribution with a half-width at 
half-maximun of about 0.6° (Figure 2, upper right panel). The 
peak of distribution is slightly shifted -0.2° from the Galactic 
plane. The distribution in Galactic longitude is noticeable based 
on the presence of many groups of clusters associated with large 
star-forming complexes (Figure 2, lower panel). For example, 
clusters located at I = 305° - 306° are associated to the star¬ 
forming complex G305; at I = 309° are associated with the Hii 
region RCW 80, clusters located at / = 310° with the Hii region 
RCW 83, at / = 326° with the star-forming complex RCW 95 
and RCW 97, at I = 333° - 334° with the star-forming com¬ 
plex RCW 106, at / = 339° with Hii region RCW 109, and the 
small maximum of cluster count located at I = 346° with the 
giant molecular complex GMC 345.5-fl.O (Lopez et al. 2011). 
We notice a decrease in the number of discovered clusters at 
I > 340°, which may be related to an increase in extinction 
and/or increased crowding toward the inner disk, which dilutes 
the ability to detect poor star clusters. 

4.2. Spatial correlation with other infrared sources and 
structures 

The methodology used to detect candidates certainly favors the 
discovery of compact groups. Detected densities were boosted 
in regions near or in the dark clouds, where the number of held 
stars drops owing to the presence of these clouds. Most of these 
cluster candidates are compact groups whose stars undergo high 
reddening, and many times they have associated small emission 
nebulosities. Some examples of this kind of clusters are La Ser¬ 
ena 045, La Serena 185, and La Serena 328. At hrst glance, the 


morphology of these star cluster candidates suggests that they 
are in the stage of ongoing star formation. Therefore, if the star 
formation processes are in such early stages, then many of these 
clusters could be associated with mid-infrared sources. 

Therefore we performed a number of cross-matching proce¬ 
dures between our catalog and the different catalogs available in 
SIMBAD and Vizier in order to identify the mid-infrared (mid- 
IR) source counterparts. The cross-matching procedures utilize 
TopCat virtual observatory tool (Taylor 2005). 
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Fig. 3: Histogram of the number of IRAS point sources associ¬ 
ated with VVV-La Serena star cluster candidates. 

The hrst two cross-matching procedures are based on the In¬ 
frared Astronomical Satellite (IRAS) point source catalog (PSC) 
(Beichman et al. 1988) and the Midcourse Space Experiment 
(MSX) catalog version 6 (MSX6C) (Egan et al. 2003). IRAS 
PSC includes about 250,000 reliable infrared point sources ob¬ 
served by the Infrared Astronomical Satellite. These sources 
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Separation [arcsec] 


Fig. 4: Histogram of the number of MSX sources associated with 
VVV - La Serena star cluster candidates. 


have angular extents less than approximately 0.5, 0.5, 1.0, and 
2.0 arcmin in the in-scan direction at 12, 25, 60, and 100yum, 
respectively. Away from confused regions of the sky, the sur¬ 
vey is complete to about 0.4, 0.5, 0.6, and 1.0 Jy at these wave¬ 
lengths. The representative position uncertainties are about 2" 
to 6" in-scan and about 8" to 16" cross-scan directions, respec¬ 
tively. The MSX6C in the Galactic plane includes about 430,000 
point sources observed with the infrared instrument SPIRIT III 
on MSX, with an astrometric accuracy better than 1". This in¬ 
strument consists of a 33 cm clear-aperture off-axis telescope 
with five line scanned infrared focal plane arrays of 18.3" square 
pixels, and a high sensitivity (0.1 Jy at 8.3 pm). SPIRIT III is pro¬ 
vided with six photometric bands: B1 (4.29yum), B2 (4.35 yum), 
A (8.28yum), C (12.13yum), D (14.65 yum), and E (21.34yum). 

We find 293 IRAS sources closer than a cross-matching ra¬ 
dius of 30" (Figure 3), corresponding to 59% of VVV cluster 
candidates. The matching radius used corresponds to the IRAS 
spatial resolution at 12yum; however, this may be too restric¬ 
tive with respect to the IRAS resolution, but will reduce the 
confusion with neighboring Mid-IR structures. The IRAS point 
sources associated with each star cluster candidate are listed in 
column 7 of Table 1. 

We also matched the cluster candidates with Midcourse 
Space Experiment (MSX) catalog version 6 (Egan et al. 2003). 
Compared with IRAS, the MSX survey has better spatial reso¬ 
lution (18") and sensitivity on the Galactic plane, and we might 
expect a tighter cross-matching. Figure 4 shows the distribution 
of separations for the MSX sources, using a matching radius of 
30". We have 436 MSX point sources (88%) associated with 
our star cluster candidates, with 50% of them inside the 8" cir¬ 
cle. The MSX sources associated with our source catalog are 
indicated in column 8 of Table 1. For sparse open clusters can¬ 
didates, composed presumably of red giant or supergiant stars 
distributed over several square arcminutes (e.g.. La Serena 314), 
the MSX point source listed in Table 1 corresponds to the closer 
source to the position derived for the cluster center. 

As mentioned above, the vast majority of cluster candidates 
were detected in areas with dark and/or bright nebulosities on the 
VVV images. A visual inspection of VVV images finds a total of 
456 cluster candidates associated with dark clouds up to a dis¬ 
tance of 2'. We searched for dark cloud catalogs in the Vizier 
database in the area covered by the Galactic disk VVV tiles, 
and we find four main entries: Peretto & Fuller (2009), Dobashi 


(2011), Dutra & Bica (2002), and Hartley et al. (1986). The 
first one is the most important for our study by far. The cross¬ 
matching between our catalog and these dark cloud catalogs 
brings 310 cluster candidates (%63) associated with 552 dark 
clouds. Column 3 in Table 2 lists the dark clouds associated with 
each cluster candidate in detail. Interestingly, most of the dark 
nebulae change their appearance and morphology with the wave¬ 
length range of the observations. As an example, we can describe 
the area surrounding the cluster La Serena 208. Figure 5 shows 
the VVV color-composite image and Spitzer/GLIMPSE 8yum 
around this cluster. In the VVV color-composite image, it is pos¬ 
sible to discern a filamentary dark nebula that causes a marked 
reddening and a decline in the star count. The Spitzer/GLIMPSE 
8 yum image at the position of La Serena 208 reveals a cluster 
of knot point sources associated with filamentary arcs that are 
deployed eastward, forming a small bubble of 1'3 diameter, cat¬ 
aloged as [CPA2006] S89 (Churchwell et al. 2006). Addition¬ 
ally, an infrared dark cloud connected with the bubble is located 
4'8 to the south. This infrared dark cloud (IRDC), known as 
SDCG32L336-0.373 (Peretto & Fuller 2009), has an angular 
size of 0'525 x 0'222. Examination of the GLIMPSE 8yum image 
suggests a connection between the IRDC and the bubble through 
a series of uncataloged dark filaments. 

In view of the importance of the rich morphological infor¬ 
mation that stands out in the Spitzer 8 pm images, we visually 
inspected the nebular structures in the GLIMPSE images at the 
positions of the star cluster candidates, which are described in 
column 11 of Table 1 The GLIMPSE survey extend approxi¬ 
mately from Galactic latitude -L2 to L2, and then only 441 
cluster candidates are projected in that area, with all of them de¬ 
tectable in some way (as point sources or extended structures) 
in the 8 pm images. Spitzer 8 pm point sources are the main fea¬ 
ture for 80 cluster candidates. Very compact nebulosity associ¬ 
ated with point sources is present in 292 cluster candidates at 
that wavelength range. These kinds of sources are described as a 
knot. In many cases (156), these point sources and/or knots are 
also associated with extended irregular nebulosities (marked as 
nebula), and there are 89 clusters where this nebulosity is bubble¬ 
shaped, which were identified as bubbles or small bubbles. Tak¬ 
ing the three different nebular emission structures cataloged into 
account, we count 361 cluster candidates with some kind of neb¬ 
ular emission associated at 8 pm, that means 82% of the sample 
in the area of GLIMPSE images. Eigure 9 depicts some exam¬ 
ples of clusters candidates as seen in GLIMPSE 8 pm images. 
The square plotted inside each image corresponds to the size of 
the VVV image charts shown in Eigure 1. 

The presence of knots and nebulosities in the Spitzer images 
suggests that the clusters could be extremely young. Taking this 
possibility into account, we explore the spatial correlation be¬ 
tween the cluster candidates and other independent tracers of 
ongoing star formation, such as masers, young stellar objects 
(YSO), extended green objects (EGO), and outflows. 

Astrophysical maser emissions are proven signatures of 
ongoing star formation, in particular methanol masers (cf. 
Ellingsen 2006; Breen et al. 2013), although we must consider 
that water or hydroxyl masers have also been detected in in star¬ 
forming regions, as well as in evolved stars or supernova rem¬ 
nants. Therefore, the presence of these maser types in the vicin¬ 
ity of a cluster could be an excellent signspot that this cluster has 
ongoing star formation. We checked the presence of masers in 
the surrounding area of candidates using the SIMBAD database 
and maser catalogs stored in Vizier. The idea of our search is 
to highlight those cluster candidates where maser emission is 
detected and then to provide a reference for the type of maser 
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Fig. 5: Area of cluster candidate La Serena 208. Left panel: 
VVV color-composite image following the same scheme as in 
Figure 1. Right panel: Spitzer/GLIMPSE 8 yum. The cluster can¬ 
didate is marked with a small red circle, and the dark cloud 
SDCG32L336 - 0.373 wth a small yellow plus sign. The field 
of view of each panel is about 9.8 x 9.0 arcmin. 


present. Table 2 lists the detected maser type, the source name 
(Column 6), and a bibliographic reference for a catalog or dis¬ 
covery paper (Column 7). 

We have found that 128 cluster candidates are associated 
with maser emissions in a radius of 60", with 50% of detec¬ 
tions (63) within the 10.5" circle (Figure 6). They correspond to 
104 methanol masers, 72 water masers, and 55 hydroxil masers. 
Detections of methanol and water masers are boosted for the effi¬ 
ciency of dedicated surveys carried out by Caswell et al. (2010), 
Caswell et al. (2011), and Green et al. (2012) for the methanol 
and Walsh et al. (2011) and Walsh et al. (2014) for water. It is 
interesting to note that 39 cluster candidates simultaneously in¬ 
clude the detection of the three type of masers (methanol, hy¬ 
droxil, and water), putting particular emphasis on these cluster 
candidates. Breen et al. (2010) point out that the simultaneous 
detection of the three type of masers perhaps indicate a special 
moment in the evolution of newly formed stars; as these authors 
mention, the methanol masers are the first species to develop 
and also the first species to disappear, whereas both water and 
hydroxil masers are more persistent species. 

YSOs are the primary tracers for star-forming regions. Cat¬ 
alogs of YSOs and YSO candidates in the area of Galactic disk 
covered by VVV survey include Mottram et al. (2007) and Yu 
& Wang (2014) (derived from MSX data), Robitaille et al. 
(2008) (derived from Spitzer/GLIMPSE data), and Chan, Hen¬ 
ning & Schreyer (1996) (derived from IRAS data). Again, we 
performed a spatial cross-matching of star cluster candidates and 
YSO candidates extracted from the mentioned catalogs, using a 
60" radius. The procedure results in the detection of 207 YSOs 
candidates in the area of 158 cluster candidates (32%), which 
are 50% of clusters with YSO candidates inside the 15" circle, 
(Figure 7). 

Other interesting infrared tracers of ongoing star formation 
are outflows, extended green objects (EGOs; Cyganowski et al. 
2008), and green fuzzies (Chambers, Jackson, Rathborne & Si¬ 
mon, R. 2009). These last two types of objects have acquired 
their names from the remarkable 4.5 jjm emission that is green 
in GLIMPSE color-composite images. EGOs and green fuzzies 
have been demostrated to be associated with young high-mass 
star forming regions (Chen et al. 2009, 2011). We searched for 
catalogs about outflows and outflows candidates in ViZier, and 
found two additional catalogs. A search for ionized outflows to¬ 
ward high-mass YSO was carried out by Guzman et al. (2012) 
through multifrequency radio contnuum observations using the 



Separation [arcsec] 


Fig. 6: Histogram of the number of masers associated with the 
new VVV - La Serena star cluster candidates. 



Separation [arcsec] 


Fig. 7: Histogram of the number of YSO candidates associated 
with the new VVV - La Serena star cluster candidates. 


Australia Telescope Compact Array, and a catalog of molecular 
outflows recompiled by Wu et al. (2004). The cross-matching 
between La Serena cluster catalog with outflows and EGOs and 
catalogs brings 73 clusters candidates associated with 90 out¬ 
flows and EGO sources inside the circle with 60" radius, where 
50% are inside 6.2" (see Figure 8). 

5. Summary 

We report the discovery of 493 new near-infrared stellar cluster 
candidates in a portion of the Galactic disk in the fourth quadrant 
using the ESO/Chile public survey VISTA variables in the Via 
Lactea. The analysis of the spatial distribution shows that the 
new clusters are concentrated near the Galactic plane and feature 
conspicuous local maxima, corresponding to large star-forming 
complexes, such as G305, RCW 95, RCW 106. Almost a half of 
the sample show a very compact appearance (250 candidates). 

The majority of the discovered cluster candidates (456, 92%) 
are associated with dark infrared nebulosities, with 310 clusters 
directly associated to 552 cataloged dark clouds. 

IRAS point sources are associated with 59% of cluster sam¬ 
ple, while for the case of MSX point sources, the cross-matching 
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Fig. 9: Examples of different cluster morphologies based on images in the 8 pm band obtained with the Infrared Array Camera 
(IRAC) on board the Spizer Space Telescope. Boxes inside each panel correspond to the cluster candidate charts showed in Figure 1. 
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Separation [arcsec] 


Fig. 8: Histogram of the number of EGO and outflow candidates 
associated with the new VVV - La Serena star cluster candi¬ 
dates. 


gives 88% of the new cluster associated to a mid-infrared source. 
In view of the importance of identiflcating mid-IR counterparts 
for the vast majority of new cluster candidates, we visually in¬ 
spected GLIMPSE 8 yum images. The morphological richness of 
the images is noteworthy. In 292 cluster candidates show knotty 
morphology at 8/im, being 361 objects associated with some 
kind of extended emission nebula at this wavelength regime. 
These morphologies suggest that most of the cluster candidates 
could be extremely young. 

This scenario of extreme youth is reinforced by the cross¬ 
matching of the new cluster candidates against masers, YSO 
candidates, outflows, and EGOs. In the first instance, 128 clus¬ 
ters candidates have associated maser emission, with 104 of 
them methanol masers, which means excellent candidates for 
ongoing massive star formation. Special mention is deserved 
by 39 cluster candidates that have associated three maser types: 
methanol, water and hydroxil. Eor YSO candidates and outflows, 
158 (32%) and 73 (15%) of new clusters candidates are closer 
than 60" from them, respectively. In particular, there is a special 
set of sixteen cluster candidates that present a clear signpost of 
star-forming activity having associated simultaneosly dark neb¬ 
ulae, young stellar objects, extended green objects, and maser 
emission. The cluster candidates are La Serena numbers: 045, 
115, 119, 161, 227, 228, 241, 283, 287, 316, 319, 391, 393, 410, 
418, and 441. 

Again, these findings suggest that a substantial proportion of 
the discovered clusters are very young and have ongoing star 
formation. This scenario is consistent with the common fea¬ 
tures of the majority of the clusters discovered independently 
by Borissova et al. (2011) and Solin et al. (2014), using the same 
VVV DRl images. 

Up to now, the vast majority of star cluster candidates dis¬ 
covered on the basis of VVV images are quite compact and gen¬ 
erally surrounded by bright or dark nebulosities. With this work, 
we think that we are getting a the deeper level of discovery us¬ 
ing the VVV survey in the Galactic disk. But much remains to 
explore in relation to the population of slightly evolved young 
clusters, where the detectable low-mass stars have left the accre¬ 
tion phase. After this phase, the clusters have dissipated much 
of the molecular cloud in which they formed, and some of the 


most massive and notable members have evolved or have been 
ejected and may show a lower concentration. These clusters are 
very difficult to separate the contaminating population of the disk 
of the galaxy. Eurther opportunities to squeeze the VVV survey 
is to exploit the stack of Ks-band images obtained during the 
current synoptic stage of the project. At present, VVV is gather¬ 
ing dozens of Ks-band images that can be the raw material for 
pursuiting further immersed and fainter stellar clusters. 
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Table 1: La Serena star cluster candidates. 


Cluster ID 

R-A./2000 

Dec/2000 

1 

b 

Tile 

IRAS source 

MSX source 

Type 

Size 

8yum morphology 

La Serena 001 

11:39:13 

-63 : 29 : 04 

294.9726 

-1.7292 

dOOl 

IRAS 

11368-6312 

G294.9725-0L7290 

c 

M 

- 

La Serena 002 

11:39:22 

-63 : 28 : 11 

294.9896 

-1.7188 

dOOl 

- 


G294.9895-01.7186 

c 

M 

- 

La Serena 003 

11:40:28 

-63 : 27 : 58 

295.1026 

-1.6779 

dOOl 

IRAS 

11380-6311 

G295.1026-01.6778 

c 

S 

- 

La Serena 004 

11:42:09 

-62 : 20 : 07 

294.9846 

-0.5365 

d039 

- 


G294.9845-00.5363 

0 

M 

P+N 

La Serena 005 

11:43:25 

-62 : 25 : 43 

295.1523 

-0.5891 

d039 

- 


G295.1523-00.5890 

0 

S 

- 

La Serena 006 

11:43:28 

-62 : 27 : 20 

295.1652 

-0.6072 

d039 

- 


G295.1651-00.6070 

0 

L 

P 

La Serena 007 

11:43:38 

-62 : 25 : 17 

295.1750 

-0.5738 

d039 

IRAS 

11412-6208 

G295.1750-00.5737 

0 

L 

K+N 

La Serena 008 

11:44:34 

-62 : 40 : 47 

295.3450 

-0.7950 

d039 

IRAS 

11421-6223 

G295.3450-00.7949 

0 

S 

sB 

La Serena 009 

11:45:04 

-63 : 17 : 44 

295.5571 

-1.3788 

dOOl 

IRAS 

11426-6301 

G295.5570-0L3787 

c 

M 

- 

La Serena 010 

11:47:05 

-62 : 38 : 31 

295.6194 

-0.6887 

d039 

- 


G295.6194-00.6886 

0 

M 

P 

La Serena Oil 

11:50:13 

-62 : 08 : 44 

295.8504 

-0.1183 

d039 

IRAS 

11478-6152 

G295.8503-00.1182 

0 

M 

P+K 

La Serena 012 

11:50:25 

-62 : 49 : 13 

296.0283 

-0.7693 

d039 

- 


- 

0 

M 

P 

La Serena 013 

11:53:12 

-62 : 30 : 18 

296.2654 

-0.3903 

d040 

IRAS 

11507-6213 

G296.2654-00.3901 

0 

M 

P+K 

La Serena 014 

11:54:15 

-62 : 51 : 41 

296.4640 

-0.7096 

d040 

IRAS 

11517-6234 

G296.4640-00.7094 

c 

M 

K 

La Serena 015 

11:55:23 

-63 : 25 : 30 

296.7131 

-1.2336 

d002 

IRAS 

11528-6308 

G296.7130-0L2335 

c 

S 

- 

La Serena 016 

11:55:53 

-62 : 54 : 35 

296.0283 

-0.7693 

d040 

- 


- 

0 

S 

P 

La Serena 017 

12:00:57 

-63 : 04 ; 04 

297.2535 

-0.7558 

d040 

IRAS 

11583-6247 

G297.2535-00.7557 

0 

s 

K+sB 

La Serena 018 

12:02:21 

-62 : 36 : 32 

297.3249 

-0.2679 

d040 

IRAS 

11597-6220 

G297.3249-00.2678 

0 

L 

P+N 

La Serena 019 

12:02:45 

-63 : 07 : 02 

297.4585 

-0.7637 

d040 

- 


G297.4585-00.7636 

0 

S 

K+N 

La Serena 020 

12:03:16 

-63 : 11 : 15 

297.5343 

-0.8266 

d040 

- 


G297.5343-00.8264 

0 

L 

K 

La Serena 021 

12:03:54 

-62 : 22 : 34 

297.4538 

-0.0094 

d040 

IRAS 

12013-6205 

G297.4537-00.0093 

c 

M 

K+P 

La Serena 022 

12:03:57 

-63 : 17 : 08 

297.6252 

-0.9044 

d040 

IRAS 

12013-6300 

G297.6252-00.9043 

c 

M 

K+N 

La Serena 023 

12:12:20 

-64 : 15 : 18 

298.6992 

-1.7068 

d003 

IRAS 

12096-6358 

G298.6992-0L7067 

0 

M 

- 

La Serena 024 

12:12:29 

-62 : 49 : 49 

298.5028 

-0.2900 

d041 

- 


G298.5028-00.2899 

0 

L 

P+N 

La Serena 025 

12:12:53 

-62 : 44 : 46 

298.5355 

-0.2051 

d041 

- 


G298.5354-00.2049 

0 

L 

P 

La Serena 026 

12:12:58 

-61 : 50 : 25 

298.4088 

+0.6938 

d079 

IRAS 

12102-6133 

G298.4087+00.6938 

0 

M 

P+N 

La Serena 027 

12:15:57 

-62 : 28 : 00 

298.8450 

+0.1221 

d079 

IRAS 

12132-6211 

G298.8449+00.1221 

0 

M 

K+N 

La Serena 028 

12:16:43 

-62 : 14 : 26 

298.8878 

+0.3599 

d079 

IRAS 

12140-6157 

G298.8877+00.3599 

0 

M 

K+P 

La Serena 029 

12:18:28 

-62 : 37 : 17 

299.1532 

+0.0086 

d042 

IRAS 

12157-6220 

G299.1531+00.0086 

c 

M 

K 

La Serena 030 

12:19:38 

-62 : 57 : 21 

299.3244 

-0.3068 

d042 

IRAS 

12169-6240 

G299.3244-00.3067 

0 

L 

K+B 

La Serena 031 

12:19:54 

-62 : 55 : 42 

299.3547 

-0.2750 

d042 

IRAS 

12171-6238 

G299.3546-00.2749 

0 

M 

K+N 

La Serena 032 

12:23:08 

-62 : 40 : 50 

299.6907 

+0.0009 

d080 

- 


G299.6906+00.0010 

0 

M 

P+B 

La Serena 033 

12:23:48 

-62 : 42 : 28 

299.7701 

-0.0061 

d080 

IRAS 

12210-6225 

G299.7700-00.0060 

c 

M 

K 

La Serena 034 

12:28:42 

-63 : 26 : 19 

300.3932 

-0.6790 

d042 

IRAS 

12259-6309 

- 

0 

M 

P 

La Serena 035 

12:28:51 

-63 : 02 : 40 

300.3770 

-0.2858 

d042 

- 


G300.3770-00.2857 

c 

M 

K+P 

La Serena 036 

12:32:22 

-62 : 41 : 20 

300.7487 

+0.1016 

d081 

IRAS 

12295-6224 

G300.7486+00.1016 

0 

M 

K+sB 

La Serena 037 

12:34:17 

-61 : 55 : 17 

300.9121 

+0.8830 

d081 

IRAS 

12314-6138 

G300.9121+00.8830 

0 

M 

K 


Notes. Only a small portion of the data is provided here, the full table is only available in electronic form at the CDS. 
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Table 2: Association 



Notes. Only a small portion of the data is provided here, the full table is only available in 


different astrophysical sources. 


EGO and outflows Masers Maser references 

CH30H [VES2000] 113651.54-631209.4 2000MNRAS.317..315V 

H20 [BSG89] 294.97-1.73 1989A&AS...77..465B 

CH30H MMB 294.990-1.719 2012MNRAS.420.3108G 


EGO G298.90-f0.36 
EGO G298.89-f0.37 


CH30H MMB 299.772-0.005 2012MNRAS.420.3108G 


ic form at the CDS. 


A&A proofs: manuscript no. Almost-500-cluster-candidates-VVV 
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Table 3: Special notes and comments for each cluster candidate. 


Cluster ID Notes 


La Serena 001 
La Serena 002 
La Serena 003 
La Serena 004 
La Serena 005 
La Serena 006 
La Serena 007 
La Serena 008 
La Serena 009 
La Serena 010 
La Serena Oil 
La Serena 012 
La Serena 013 
La Serena 014 
La Serena 015 
La Serena 016 
La Serena 017 
La Serena 018 
La Serena 019 
La Serena 020 
La Serena 021 
La Serena 022 
La Serena 023 
La Serena 024 
La Serena 025 
La Serena 026 
La Serena 027 
La Serena 028 
La Serena 029 
La Serena 030 
La Serena 031 
La Serena 032 
La Serena 033 
La Serena 034 
La Serena 035 
La Serena 036 
La Serena 037 


1.6’ from cluster VVV-CL 005. 1.8’ from IC2944 and 1’ from IC2948 
1.4’ from cluster VVV-CL 101 
4’ from IC 2944 

In HII region [WMG70] 295.2-00.6 
Projected on open cluster Stock 14 

5’ from open cluster Stock 14, and 1.3’ from SNR G295.2-00.6 


Identified also as cluster FSR 1595 


HII region GAL 297.51-00.77, 3’ from cluster Mercer 29, in the edge of bubble 
3.5’ from cluster Mercer 29 

In molecular cloud [RC2004] G297.6-0.9-f 30.3 

In the edge of large bubble, cluster Mercer 30 inside. 

2.6’ from cluster DSH J1213.0-6242 

HII region. It is not PN G298.4-1-00.6. Refiection nebula GN 12.10.3 


Dark nebula SDC G299.167-f 0.010 

Refiection nebula BRAN 386C, in the wall of bubble [CPA2006] SI75. Part of RCW 64. 
3.4’ from cluster VVV-CL 012. Part of RCW 64. HII region BRAN 386E 
In the center of bubble [CPA2006] S174 

In the edge of bubble [CPA2006] S174, 4.8’ from cluster VVV-CL 032 


Notes. Only a small portion of the data is provided here, the full table is only available in electronic form at the CDS. 
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